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HIGHLIGHTS 


•  Severe  non-uniformities  in  a  DFMC  are  predicted  under  ultra-low  fuel  stoichiometry  condition. 

•  This  paper  identifies  the  controlling  parameters  of  the  methanol  transport  in  a  DMFC. 

•  Two  strategies  called  ‘streamline  graded  structures  (SGSs)’  to  mitigate  non-uniformity  are  proposed. 

•  Computational  results  show  SGS  models  have  10%  improved  performance  compared  to  the  baseline  model. 
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Operating  a  direct  methanol  fuel  cell  (DMFC)  with  ultra-low  fuel  stoichiometry  is  desirable  to  build  a 
compact  portable  power  system.  Through  3D  computational  modeling  with  statistical  analysis,  severe 
non-uniformity  of  fuel  concentration  leading  to  high  methanol  crossover  in  the  inlet  region  and  poor  fuel 
supply  in  the  outlet  region  is  predicted  under  ultra-low  fuel  stoichiometry  condition,  which  causes  large 
voltage  loss  or  even  cell  shut-down.  After  identifying  controlling  parameters  of  methanol  transport,  this 
paper  proposes  streamline-graded  structures  (SGS)  of  the  anode  channel  to  mitigate  anode  fuel  distri¬ 
bution  non-uniformity  and  to  boost  cell  performance  and  fuel  efficiency  together.  Computational  results 
show  that  streamline-graded  structures  achieve  about  10%  of  voltage  gain  and  3%  of  fuel  efficiency 
improvement  compared  to  the  conventional  DMFC  design  by  mitigating  anode  non-uniformity. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  order  to  develop  a  compact,  portable  direct  methanol  fuel  cell 
(DMFC)  system,  there  are  several  objectives  to  achieve,  such  as  high 
energy  density,  small  system  size  and  high  fuel  efficiency.  First, 
energy  density  of  the  DMFC  system  must  be  increased  by  using 
highly  concentrated  fuel.  Many  DMFCs  carry  a  diluted  fuel  solution 
which  consists  of  water  and  methanol.  Although  cell  operation  with 
dilute  solution  has  benefits  of  reducing  methanol  crossover  and 
increasing  cell  voltage,  it  is  not  optimal  for  development  of  portable 
DMFCs  because  dilute  solution  lowers  energy  density  of  the  entire 
DMFC  system.  As  water  is  produced  in  the  cathode  catalyst  layer  by 
oxygen  reduction  reaction  (ORR),  it  is  theoretically  possible  to 
achieve  a  water-neutral  condition  to  operate  a  DMFC  without 
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external  water  supply  to  the  anode  side.  Blum  et  al.  [1]  first  pro¬ 
posed  this  water-neutral  condition  and  they  concluded  that  car¬ 
rying  pure  liquid  methanol  is  theoretically  possible  when  the  water 
transfer  coefficient  (a)  is  equal  to  -1/6.  However,  supplying  pure 
methanol  to  a  DMFC  causes  severe  methanol  crossover  and  it  is 
difficult  to  uniformly  supply  fuel  to  a  DMFC  with  low  fuel  stoichi¬ 
ometry.  Even  though  a  DMFC  system  carries  a  pure  methanol  tank, 
it  must  dilute  pure  methanol  internally  with  produced  water. 

Balance  of  plant  (BOP),  occupying  a  large  volume  of  a  DMFC 
system,  should  be  reduced  also.  The  gas-liquid  separator  located  at 
the  anode  outlet  can  be  removed  if  CO2  gas  is  vented  through 
alternate  routes.  Conventional  DMFCs  use  a  pump  to  supply  liquid 
fuel  on  the  anode  side.  If  anode  flow  rate  is  minimized,  this  pump 
size  can  be  reduced  and  required  power  for  pump  operation  can  be 
reduced  as  well.  Due  to  C02  production  by  methanol  oxidation 
reaction  (MOR)  in  the  anode,  a  strong  two-phase  flow  exists  in  the 
anode  channel.  This  two-phase  flow  causes  a  large  pressure  drop  by 
two-phase  frictional  loss  and  inertial  loss.  The  C02  breathing  type 
DMFC  removes  C02  directly  to  the  ambient  and  almost  pure  liquid 
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flows  in  the  anode  channel,  which  reduces  the  required  pumping 
power  [2,3]. 

Maintaining  ultra-low  fuel  stoichiometry  is  another  way  to 
construct  a  compact  DMFC  system  and  reduce  manufacturing  cost. 
When  a  DMFC  operates  with  large  fuel  stoichiometry  condition, 
redundant  methanol  which  did  not  participate  in  anode  reaction  is 
lost  to  the  cathode  side  by  diffusion  and  electro-osmosis  drag 
through  the  proton  exchange  membrane.  This  methanol  crossover 
lowers  the  overall  fuel  efficiency.  High  fuel  efficiency  is  particularly 
important  for  a  portable  system  as  it  determines  the  total  opera¬ 
tional  time  of  the  system.  The  easiest  way  to  minimize  methanol 
crossover  is  to  operate  a  DMFC  with  ultra-low  fuel  stoichiometry. 


V-(ypt?Yk)  =  V-  j/>,DkeffVY|k  +  pgDgeffVYk] 

-v[(''li-',i)T,]+sk  (,) 

When  the  vapor-liquid  equilibrium  of  methanol  is  assumed  in 
the  anode  by  Flenry’s  law,  gas  and  liquid  methanol  concentration 
are  related  as 

cMeOH  =  cMeOH^MeOH 

The  following  relationship  can  be  established  by  multi-phase 
mixture  model  theory. 


pyMeOH 


piyMeOHs  _j_  ^yMeOH  -  s)  where 


p  yMeOH  =  cMeOHj^MeOH 
p^yMeOH  =  cMeOHj^MeOH 
pgyMe°H  =  (!MeOHMMeOH 


(3) 


And,  theoretically,  it  is  ideal  to  operates  a  DMFC  near  anode  limiting 
current  density  regime  where  methanol  crossover  is  almost  zero, 
which  enables  the  cell  to  achieve  maximum  performance  with 
maximum  fuel  efficiency.  It  is,  however,  not  easy  to  achieve  this 
ideal  operating  condition  because  fuel  concentration  may 
distribute  in  three  dimensions  non-uniformly. 

When  a  DMFC  operates  at  very  low  fuel  flow  rate,  convectional 
effect  of  methanol  transport  in  the  channel  direction  becomes 
weak.  Instead,  diffusive  effect  of  methanol  transport  through  the 
membrane  electrode  assembly  becomes  stronger  and  large  amount 
of  methanol  may  diffuse  through  the  proton  exchange  membrane 
in  the  inlet  region,  causing  severe  methanol  crossover  there.  In 
contrast,  the  outlet  region  may  suffer  from  fuel  depletion,  which 
also  degrades  cell  performance  by  increasing  mass  transport 
overpotential.  Therefore,  when  liquid  fuel  is  supplied  from  the  inlet 
at  ultra-low  flow  rate,  some  regions  will  operate  with  excessive  fuel 
and  some  regions  may  suffer  from  fuel  shortage.  This  is  called  ‘fuel 
concentration  non-uniformity’.  Fuel  concentration  non-uniformity 
eventually  leads  to  current  density  non-uniformity,  temperature 
non-uniformity,  and  even  cell  shut-down.  All  of  those  non¬ 
uniformities  are  called  ‘anode  non-uniformity’.  It  can  be  easily 
conjectured  that  the  anode  non-uniformity  has  a  strong  influence 
on  the  overall  cell  performance  and  fuel  efficiency. 

In  this  paper,  we  mathematically  investigate  the  anode  non¬ 
uniformity  as  well  as  its  consequences  on  cell  performance  when 
a  DMFC  operates  at  ultra-low  fuel  stoichiometry.  After  identifying 
that,  strategies  of  mitigating  anode  non-uniformity  are  introduced 
in  order  to  improve  both  cell  performance  and  fuel  efficiency 
required  for  portable  power  sources. 


2.  Physics  model 

In  order  to  analyze  transport  phenomena  in  a  DMFC,  previously 
developed  multi-D  DMFC  model  [4—6]  is  applied  to  the  present 
study.  Governing  equations  and  constitutive  relationships  are 
summarized  in  Tables  1  and  2,  respectively.  Coupled  PDEs  (mass, 
momentum,  species,  proton  transport,  electron  transport,  and  en¬ 
ergy)  are  solved  together  with  source  terms.  Cell  dimension  and 
simulation  parameters  are  presented  in  Tables  3  and  4,  respectively. 
As  the  present  study  focuses  on  methanol  transport,  the  methanol 
transport  equation  for  the  multi-D  model  is  revisited  here.  A  gen¬ 
eral  form  of  multi-phase  species  equation  by  M2  model  [7]  can  be 
expressed  as  follows: 


Combining  Eq.  (2)  and  Eq.  (3),  methanol  mass  fraction  can  be 
found  as  follows 


yMeOH 


A 

P 


1  -S 
kfp H 


yMeOH 


(4) 


Therefore,  the  methanol  transport  equation  can  be  reformed 
from  Eq.  (1)  as  follows: 


p-yMe°H) 


=  V- 
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rjMeOHx 

nMeOH  ,  g,eff  \  ^ 
U\,eff  +  j^MeOH  J  V 


I  \7  •  VMe0H 


-  V 
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yMeOH 


-Y 


MeOHV 
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(5) 
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MeOH 


The  first  term  on  the  right  hand  side  in  Eq.  (5)  expresses  the 
two-phase  molecular  diffusion.  Since  molecular  diffusion  of 
methanol  is  carried  out  by  both  liquid  phase  and  gas  phase,  we 
define  two-phase  molecular  diffusivity  of  methanol  as  follows: 


D 


|MeOH 

eff 


=  ^  D 


'jMeOH 

3,eff 
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DMeOH 

g,eff 

^MeOH 


=  —\  D\ 


'»MeOHcn 


DMe°H(i  _S)H 


kMeOH 


(6) 


As  the  value  of  /<^e0H  is  about  1191  in  the  room  temperature, 
Dj^pH  and  D^|^H//<ffe0H  becomes  similar  order  of  magnitude.  Note 
that  two-phase  molecular  diffusivity  of  methanol  is  a  strong  func¬ 
tion  of  both  operating  condition  (liquid  saturation)  and  material 
property  (porosity).  The  liquid  saturation  shown  above  equation  is 
defined  as  follows  [7]: 


.  vi  1  Vg  pgyH2°-pgygH2° 

Vpore  Vpore  PxY^°  -  pxYH^°  +  pgY^O  _  p^O 

(7) 


And  porosity  is  defined  as, 


=  Ypore  ^pore 

V  Vg  +  Vj  +  1/pore 

The  second  term  in  Eq.  (5)  explains  methanol  transport  by 
capillary  force  which  becomes  significant  when  the  anode  flow  rate 
is  large  or  methanol  concentration  is  high.  As  the  present  study 
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Table  1 

Governing  equations  of  the  multi-D  DMFC  model. 


Conservation  equations 

Source  terms 

Mass 

V-(pT?)  =  Sm 

S™d  =  MMe0HSMe0H  +MH2°SH2°  +  ^Mc°2 

S™d  =  mH2°SH2°  +M°2S°2  +ipiMc°2 

Momentum 

v-(pTtit)  =  -vp  +  v-V  +  Su 

Su  =  -j (in  porous  media  domain) 

Water  transport 

v.(T/,lfyH20)  =  v-[/>DcapUIVYH20]+SH!° 

s£°  =  -^MH2°  -JTaMH2° 

S^o  =  ^MH2°  +%tMH2°+f«MH2° 

Methanol  transport 

V- (-)-/>  U>Y|Me0H)  =  V  ■  [/)D£{f0H  V  •  yMeOH]  +  sMe0H 

^MeOH  _  j  jy^MeOH  _  jxoverjy^MeOH 

Oxygen  transport 

V(jptY^)  =  V-^D^VY02]  +V-[Y°27,]  +S°2 

s°2i  =  “iM°2  -VM°2 

Proton 

0  =  v-(Ke>effv^e)  +  s* 

5 0  —  j 

^acl  ~  J 

•^e.ccl  =  —jc  +fx  over 

Electron 

o  =  v'(Ks,effV0s)  +sf 

5 0  —  —i 

s,acl  “  J 

stcd  =  jc  —jx  over 

Energy 

V(pcpuT)  =  V(keffVT)+ST 

Sensible, acl  =j(v  +  T^jf)  +  ^ 

cT  —  JL 

Jsensible,mem  KeS 

^sensible, ccl  =  j(v  +  7 IFf)  +  iw  +  Jxovel  (v  +  T +) 
SLent  =  V-(h°  ftT?,) 

focuses  on  ultra-low  stoichiometry  condition,  the  anode  becomes 
gaseous  and  the  methanol  capillary  term  can  be  neglected.  Finally, 
the  methanol  transport  equation  becomes, 

V  •  (Vy,Me0H)  =  V  •  [pD“e0H  V  •  y,Me0H]  +  sMe0H  (9) 

Fig.  1  schematically  shows  methanol  transport  in  the  anode. 
Supplied  methanol  from  the  channel  inlet  spreads  out  through  the 


anode  porous  media  and  the  remaining  amount  moves  out  of  the 
channel  outlet.  Methanol  is  convection-dominantly  transported 
from  the  channel  to  the  surface  of  the  gas  diffusion  layer.  Then, 
methanol  is  diffusively  transported  to  the  catalyst  layer  through  the 
porous  media.  Note  that  the  ratio  of  channel  width  and  membrane 
width  strongly  affects  diffusive  mass  transport  between  the  surface 
of  the  gas  diffusion  layer  and  the  catalyst  layer.  Considering  the 
above  two  mass  transport  mechanisms,  methanol  transported  from 
the  channel  to  the  catalyst  layer  can  be  simply  expressed  as  Eq.  (10). 


Table  2 

Constitutive  relationships  and  parameters  for  the  multi-D  DMFC  model. 


Parameters 

Expression 

Reference 

Relative  permeability 

fcr.l  =s":  kr,g=(1-s)n 

[14] 

Methanol  diffusivity  (vapor) 

DMeOH,H20  =  DMeOH,C02  =  j  96  x  10-5  (^r^'  ^loUOO.O 

[15] 

Water  diffusivity  (vapor) 

DH20,C02  =  2  Q1  x  1 0“5  (3J?) 1  823 101300.0 

[16] 

Methanol  diffusivity  (liquid) 

DMe0H,H20  =  1  4  x  10-9(6473-298.15)<i 

[20] 

Oxygen  diffusivity  (gas) 

Dg2  =  3.57  x  10-5(332)’  823,013000 

[16] 

Water  content  in  Nafion-membrane 

7  -  |  22  !  gffjf  (S  <  0.3)0.043  +  17.81 RH  -  39.85 RH2  +  36.0 RH3  (vapor) 

[17,18] 

Water  diffusivity  in  Nafion-membrane 

Dft 'em  =  4.80  x  10-"exp[2416(3J3- ')] 

Calibrated  based  on  [17] 

Methanol  diffusivity  in  Nafion-membrane 

«  =  1.5xlO-,0exp[2416(+-})] 

Calibrated  based  on  [17] 

Henry’s  constant 

~  101.235  exP[4210(r  298.15)]  r 

Fitted  data 

Water  EOD  coefficient  in  Nafion-membrane 

nH2°  =  |  (nJ2°f  _  -L)  +  !  (for  X  >  14)1.0  (for  A  <  14) 

[19],  [21] 

Water  reference  EOD  coefficient  in  Nafion-membrane 

nd2r?f  =  1-6767 +  0.0155(1 -273)  + 8.9074  X  10“5(T-273)2 

[19] 

Methanol  drag  coefficient  in  Nafion-membrane 

„MeOH  _  nH2OcMeOH 
nd  -  nd  “c^T 

[5] 

Ion  conductivity  of  Nafion-membrane 

/<mem=0.1(Scm_1)@60°C  (fully  hydrated) 

Effective  water  vapor  diffusivity 

Dg.eff  =  Ds2°(1  -s)"£n 

J-Leverett  function 

,/  x  f -1.263s3  +  1.669s2  -0.966s +  0.56  (when  dc  <  90°) 

J[  )  \  1.263s3  -  2.120s2  +  1.417s  (when  dc  >  90°) 
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Fig.  1.  Schematic  of  methanol  transport  by  convection  and  diffusion  in  the  anode  of  a 
DMFC. 


•MeOH  _  Anem  { ^ea°H  C^0H  }  MMeQH 
_ 1 _ | _ _ 

1 1 _2— 1  /n ^ -Ar) j-  (10) 

where  Ar  =  (openning  ratio) 

Note  DMe°H{i  —  2~1/n(l  —  Ar)}  in  the  above  equation  is  a  semi- 
empirical  equation  derived  for  considering  mass  transport  from 
the  surface  of  the  anode  gas  diffusion  layer  to  the  catalyst  layer. 

The  effective  resistance  which  consists  of  convective  resistance 
and  diffusive  resistance  expresses  the  total  resistance  controlling 
methanol  transport  as  follows: 

Reff  =  Rconv  +Rdiff  =  Sn^  +  DMeOH{1_2-Vn(l-/lr)} 

(ID 

Convective  resistance  (RConv)  is  determined  by  opening  ratio  (Ar) 
and  convective  mass  transfer  coefficient  which  is  a  function  of  flow 
rate  [Re d)  as  follows  [8]: 


as'o,aC“fOHexp(^J?a) 
cMe°H  +  kaexp(^7Ja) 


where  rja  =  0S  -  0e  -  U ° 


(13) 


Eq.  (13)  explains  the  reaction  order  shift  from  1st  to  0th  when 
methanol  concentration  becomes  higher  than  the  threshold  value. 
On  the  cathode  side,  oxygen  reduction  reaction  (ORR)  is  assumed  to 
be  proportional  to  oxygen  concentration. 


jc  =  as<o,cccd(l  -s)nexp  where  rjc 

=  0s-0e-U°  (14) 

Fig.  3(a)  shows  the  base  model  geometry  used  in  the  present 
study.  Micro-porous  layers  are  inserted  between  gas  diffusion  layer 
and  catalyst  layer  on  both  anode  and  cathode  side.  It  is  known  that 
saturation  jump  occurs  between  different  porous  media,  and  the 
micro-porous  layer  plays  an  important  role  in  controlling  water 
management  in  a  DMFC  by  developing  saturation  jump  [12,13]. 
Previous  study  [6]  mathematically  revealed  that  the  anode  micro- 
porous  layer  acts  as  a  barrier  which  preventing  liquid  in  the 
anode  from  flowing  to  the  cathode,  which  helps  to  lower  methanol 
crossover  and  water  transfer  coefficient,  whereas  the  cathode 
micro-porous  layer  reserves  liquid  in  the  cathode  catalyst  layer, 
increasing  water  back-diffusion.  Multi-D  saturation  jump  model  is 
applied  to  the  present  DMFC  model  which  has  four  interfaces  be¬ 
tween  different  porous  media.  The  anode  becomes  gaseous  in  the 
channel  direction  since  CO2  is  accumulated,  whereas  the  cathode 
becomes  wetter  in  the  channel  direction  since  produced  water  is 
accumulated  according  to  the  operational  condition  (cathode  inlet 
gas  humidity).  In  addition,  the  region  under  the  anode  land  is 
gaseous  due  to  CO2  removal  blockage  there,  whereas  the  region 
under  the  cathode  land  is  wet  due  to  water  removal  blockage.  The 
energy  equation  is  solved  in  order  to  examine  temperature  distri¬ 
bution  which  may  be  affected  by  anode  non-uniformity  in  the 
present  study. 

The  present  multi-D  model  is  implemented  into  a  commercial 
CFD  package,  STAR-CD™  through  the  user  subroutine  capability. 
AMG  (algebraic  multi-grid  method)  based  on  finite  volume  method 
(FVM)  is  used  for  discretization  and  SIMPLE  algorithm  is  applied  for 
solving  governing  equations.  It  is  considered  that  convergence  is 
achieved  when  the  residuals  reach  10-8. 


h  m 


Sh 


DcMhea°nH 


H 


1.86 


( R^ScN 
1  L-Dh  ) 


1/3 


pan 


where 


l  ReD„ 


PuDh 

P 


(12) 


Under  ultra-low  fuel  stoichiometry  condition,  methanol 
diffusive  resistance,  Rdiff  becomes  the  dominant  transport  resis¬ 
tance,  which  is  a  function  of  both  material  property  (porosity) 
and  channel  geometry  (Ar).  Dependence  of  the  methanol 
diffusive  resistance  on  porosity,  channel  geometry  and  temper¬ 
ature  is  presented  in  Fig.  2.  Diffusive  resistance  exponentially 
increases  when  porosity  decreases.  Opening  ratio  and  tempera¬ 
ture  also  affect  the  diffusive  resistance,  albeit  less  strongly  than 
porosity. 

The  effect  of  methanol  concentration  in  the  anode  catalyst  layer 
on  the  cell  performance  is  considered  as  follows  [9—11]: 


3.  Results  and  discussion 

3.1.  Anode  fuel  non-uniformity  under  ultra-low  flow  stoichiometry 

The  effect  of  anode  stoichiometry  on  the  cell  performance  and 
fuel  non-uniformity  is  investigated  using  the  multi-D  simulation. 
Methanol  concentration  distribution  in  the  anode  catalyst  layer  for 
stoichiometry  2.4  and  1.4  are  presented  in  Fig.  4(a)  and  4(b), 
respectively.  Under  ultra-low  anode  stoichiometry  condition 
(?fl  =  1.4),  methanol  is  concentrated  in  the  inlet  region  and  the 
outlet  half  region  suffers  from  fuel  shortage.  Current  density  or 
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Fig.  2.  Diffusive  resistance  according  to  (a)  GDL  porosity  (c),  (b)  opening  ratio  (Ar),  and  (c)  Temperature  (K). 


anode  reaction  rate  is  strongly  affected  by  methanol  concentration 
when  methanol  concentration  is  insufficient.  For  the  same  oper¬ 
ating  current  condition  (150  mA  cm-2),  the  outlet  half  region  has 
lower  current  density  than  the  average  value  whereas  the  inlet  half 
region  has  high  current  density  compared  to  the  average  value, 
which  means  large  anode  non-uniformity  of  current  density  (see 
Fig.  5).  Crossover  current  density  which  is  calculated  by  Eq.  (15)  is 
presented  in  Fig.  6. 

Ixover  P= 

whorp  /VMe0H  -  0MeQHC^'°H  |  „MeOHi 
wneie  ivmem  —  umem  -h  nd  f 

As  methanol  crossover  consists  of  diffusion  through  the  mem¬ 
brane  and  electro-osmosis  by  electric  current,  methanol  crossover 
is  concentrated  in  the  inlet  region  where  both  current  density  and 
methanol  concentration  are  high  under  ultra-low  anode  stoichi¬ 
ometry  condition.  When  fuel  stoichiometry  is  large,  methanol 
crossover  becomes  more  uniform.  Methanol  crossover  rate  also 


affects  temperature  distribution.  Examining  heat  source  terms  in 
the  energy  equation  (see  Table  1),  crossover  current  density  con¬ 
tributes  to  the  heat  generation  on  the  cathode  side.  As  methanol 
crossover  is  concentrated  in  the  inlet  region,  the  low  anode  stoi¬ 
chiometry  case  shows  larger  temperature  gradient  along  the 
streamline  direction  compared  to  the  large  stoichiometry  case  (see 
Fig.  7(a)  and  7(b)). 

Although  large  anode  flow  stoichiometry  assures  better  cell 
voltage  due  to  lowered  mass  transport  overpotential  and  good 
anode  uniformity,  it  causes  large  methanol  crossover  rate  and  re¬ 
quires  considerable  fuel  pumping  power.  In  contrast,  the  cell  may 
suffer  from  severe  fuel  shortage  at  the  outlet  region  under  ultra- 
low  anode  stoichiometry  condition,  which  results  in  large  anode 
overpotential  and  poor  anode  uniformity. 

High  temperature  leads  to  better  methanol  diffusion  in  porous 
media  due  to  the  Arrhenius  relationship  of  diffusivity  as  shown  in 
Fig.  2(c).  In  addition,  methanol  oxidation  reaction  is  also  improved 
from  Eq.  (13).  Therefore,  high  temperature  generally  leads  to 
better  cell  performance.  However,  this  may  not  be  the  case  when 
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Fig.  3.  Three-dimensional  geometry  of  DMFC  models  which  have  seven  porous  layers 
(anode  diffusion  layer  |  anode  micro-porous  layer  |  anode  catalyst  layer  |  proton  ex¬ 
change  membrane  |  cathode  catalyst  layer  |  cathode  micro-porous  layer  |  cathode 
diffusion  layer),  (a)  Base  DMFC  (upper  model),  (b)  SGS  DMFC  with  tapered  channel 
(lower  model). 


the  cell  operates  at  ultra-low  fuel  stoichiometry.  High  temperature 
increases  not  only  the  effective  diffusivity  in  the  anode  porous 
media  but  also  methanol  diffusivity  of  the  proton  exchange 
membrane,  which  is  not  desired  for  fuel  management  since  it  may 
exacerbate  large  methanol  crossover.  Effect  of  temperature  on  fuel 
distribution  can  be  examined  by  comparing  Fig.  4(b)  with  (c).  As  a 


large  amount  of  methanol  is  transported  to  the  anode  catalyst 
layer  in  the  inlet  region  due  to  increased  effective  methanol 
diffusivity,  much  more  methanol  is  consumed  there  by  both 
methanol  oxidation  reaction  and  methanol  crossover  when  the 
temperature  is  high  under  ultra-low  anode  stoichiometry  condi¬ 
tion.  Therefore,  the  outlet  half  region  suffers  from  fuel  shortage 
and  poor  anode  kinetics,  which  worsens  anode  non-uniformity 
(see  Figs.  5-7).  It  is  predicted  that  the  cell  may  shut  down  with 
severe  anode  non-uniformity  when  the  cell  temperature  exceeds 
70  °C  under  ultra-low  anode  stoichiometry  in  the  present  study. 
Therefore,  excessive  cell  temperature  should  be  managed  under 
ultra-low  anode  stoichiometry  condition.  In  order  to  operate  the 
cell  at  high  temperature  while  maintaining  low  methanol  cross¬ 
over  under  ultra-low  anode  flow  stoichiometry  condition,  a  hy¬ 
drocarbon  membrane  with  low  methanol  diffusivity  can  be 
selected  as  the  proton  exchange  membrane  instead  of  a  fluoro¬ 
carbon  membrane. 


3.2.  Mitigation  of  non-uniformity  by  controlling  methanol 
transport  resistance 

It  is  clear  that  non-uniform  fuel  distribution  leads  to  other 
anode  non-uniformities  such  as  current  density,  crossover  current 
density,  and  temperature  as  discussed.  In  addition,  severe  anode 
non-uniformity  leads  to  poor  cell  voltage  and  poor  fuel  efficiency. 
Therefore,  if  we  mitigate  methanol  concentration  non-uniformity, 
we  can  simultaneously  improve  cell  performance  and  fuel  effi¬ 
ciency.  Cell  operation  under  ultra-low  anode  stoichiometry  can 
thus  be  achieved  without  sacrificing  cell  performance.  In  order  to 
mitigate  non-uniform  fuel  distribution,  we  should  investigate 
clues  for  this  in  the  methanol  transport  mechanism  discussed 
previously. 


Fig.  4.  Methanol  concentration  distribution  in  the  anode  catalyst  layer  of  DMFCs  operating  at  i  =  150  mA  cm  2,  (a)  Base  model  (Tceii  =  313  K,  £a  =  2.4,  £c  =  2.4),  (b)  Base  model 
(TCeii  =  313  K,  Ca  =  1-4,  Cc  =  2.0),  (c)  Base  model  (Tceii  =  323  K,  £a  =  1.4,  £c  =  2.0),  (d)  SGS  model  with  variable  porosity  type  A  (Tceii  =  313  K,  £a  =  1.4,  £c  =  2.0),  (e)  SGS  model  with 
tapered  channel  (Tceii  =  313  K,  £a  =  1.4,  £c  =  2.0). 
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Fig.  5.  Current  density  distribution  of  DMFCs  operating  at  i  =  150  mA  cm-2,  (a)  Basel  model  (Tceii  =  313  K,  Ca  =  2.4,  £c  =  2.4),  (b)  Base  model  (Tceii  =  313  K,  £a  =  1.4,  £c  =  2.0),  (c)  Base 
model  (TCeii  =  323  K,  Ca  =  1-4,  £c  =  2.0),  (d)  SGS  model  with  variable  porosity  type  A  (Tcen  =  313  K,  £a  =  1.4,  £c  -  2.0),  (e)  SGS  model  with  tapered  channel  (Tcen  =  313  K,  £a  =  1.4, 
Cc  -  2.0). 


It  is  shown  that  the  effective  resistance  in  Eq.  (11)  controls 
methanol  transport  in  a  DMFC.  The  conventional  DMFC  has 
almost  uniform  effective  methanol  transport  resistance  in  the 
streamline  direction,  which  results  in  large  methanol  concen¬ 
tration  gradient  in  the  streamline  direction,  i.e.,  methanol  con¬ 
centration  non-uniformity  (see  Figs.  4-7).  In  order  to  build  a 
more  uniform  methanol  concentration  distribution  in  the 
streamline  direction,  we  should  vary  the  effective  methanol 
transport  resistance  itself.  The  inlet  region  where  large  methanol 
crossover  occurs  must  have  significant  methanol  transport 
resistance  to  reduce  methanol  transport.  In  contrast,  the  outlet 
region  must  have  minimal  methanol  transport  resistance  to 
enhance  methanol  transport.  Hereafter,  we  call  this  concept  the 
streamline-graded  structure  (SGS).  Although  cell  temperature 
and  methanol  concentration  also  affect  the  effective  methanol 
diffusivity,  they  are  not  design  parameters  but  operating  condi¬ 
tions.  Opening  ratio  and  porosity  are  important  design  parame¬ 
ters  to  construct  SGS,  and  affect  the  effective  methanol  transport 
resistance  as  shown  in  Fig.  2. 

The  opening  ratio  affects  not  only  convective  transport  but  also 
diffusive  transport.  Narrow  channel  width  (small  Ar)  in  the  inlet 
region  reduces  diffusive  transport  by  increasing  diffusion  length. 
Deepening  channel  depth  can  be  helpful  for  reducing  convective 
effect  and  delaying  methanol  diffusion  in  inlet  region.  On  the 
contrary,  in  the  outlet  half  region  where  methanol  concentration  is 
low,  a  wide  channel  is  beneficial  due  to  short  diffusion  length,  and 
shallow  depth  improves  the  convective  effect.  In  summary,  the  final 
anode  channel  design  has  a  three-dimensionally  tapered  shape  as 
shown  in  Fig.  3(b). 


Molecular  methanol  diffusivity  (D^^H)  is  a  strong  function  of 
porosity.  In  order  to  get  a  linearly  varying  diffusivity  in  streamline 
direction  as, 


nMeOH{.,\  _  nMeOH 
^mol  O')  -  ^min 


D^OH_DMeOHX 
max _ min  j . , 

1  J 


(16) 


Porosity  of  the  anode  gas  diffusion  layer  can  be  set  up  as, 


(17) 


However,  it  is  not  easy  to  fabricate  a  material  which  has 
smoothly  varying  porosity  as  above.  Instead,  porosity  can  be 
distributed  in  a  stepping  shape  as  shown  in  Fig.  8. 


3.3.  Methanol  transfer  resistance 

Liquid  saturation  in  the  anode  of  the  multi-D  model  decreases 
along  the  channel  due  to  CO2  gas  accumulation.  Therefore,  mo¬ 
lecular  methanol  diffusivity,  which  is  a  strong  function  of  liquid 
saturation,  increases  and  methanol  transport  resistance  decreases 
in  the  channel  direction  in  the  base  model  as  shown  in  Fig.  9.  SGS 
models  show  much  greater  resistance  in  the  inlet  region  compared 
to  the  base  model.  Although  each  SGS  model  adopts  different 
approach  to  control  methanol  transport,  the  resulting  resistance  is 
similar. 
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Fig.  6.  Crossover  current  density  distribution  of  DMFCs  operating  at  i  =  150  mA  cm-2,  (a)  Base  model  (Tceii  =  313  K,  £a  =  2.4,  £c  =  2.4),  (b)  Base  model  (Tceii  =  313  K,  £a  =  1.4, 
Cc  =  2.0),  (c)  Base  model  (Tceu  =  323  K,  Ca  =  1.4,  Cc  =  2.0),  (d)  SGS  model  with  variable  porosity  type  A  (Tceii  =  313  K,  £a  =  1.4,  Cc  =  2.0),  (e)  SGS  model  with  tapered  channel 
(Teen  =  313  K,  Ca  -  1-4,  Cc  =  2.0). 


3.4.  Methanol  concentration  distribution 

It  is  important  to  avoid  1  st-order  methanol  oxidation  reaction  in 
the  anode  catalyst  layer  for  high  uniformity  of  current  density 
while  maintaining  sufficiently  low  methanol  concentration  near 
the  anode  catalyst  layer  to  avoid  severe  methanol  crossover  which 
deteriorates  both  cell  performance  and  fuel  efficiency.  This 
requirement  eventually  leads  to  achieving  high  uniformity  of 
methanol  concentration  distribution,  which  is  the  main  goal  of  the 
present  study. 

Methanol  concentration  in  the  anode  channel  direction  is  pre¬ 
sented  in  Fig.  10.  Contrary  to  the  base  model  which  shows  large 
variation  along  the  channel,  SGS  models  show  relatively  uniform 
distribution  of  methanol  concentration  due  to  the  effect  of  redis¬ 
tribution  of  methanol  transport  resistance.  Compared  to  the  base 
model  (Fig.  4(b)),  SGS  models  (Fig.  4(d)  and  (e))  shows  increased 
methanol  concentration  in  the  outlet  region  while  it  is  lowered  in 
the  inlet  region.  This  mitigates  anode  non-uniformity  that  will 
reduce  methanol  crossover  in  the  inlet  region  and  risk  of  fuel 
shortage  in  the  outlet  region  also. 

3.5.  Current  density  distribution 

Calculation  results  of  current  density  distribution  in  the  channel 
direction  of  multi-D  model  are  presented  in  Fig.  11.  As  the  cathode 
side  is  same  for  all  cases,  different  current  density  distribution 
comes  from  the  anode  side.  The  base  model  shows  severe  non¬ 
uniformity  which  is  direct  consequence  of  non-uniform  distribu¬ 
tion  of  methanol  concentration. 


In  the  base  model,  peak  current  occurs  not  just  in  the  inlet 
region  but  around  the  center  region  in  the  channel  direction  due 
to  large  methanol  crossover  at  the  inlet.  SGS  models  have 
improved  uniformity,  and  peak  current  density  is  observed 
almost  in  the  inlet  region  due  to  reduced  methanol  crossover. 
A  contour  plot  of  current  density  is  shown  in  Fig.  5.  Note  that 
peak  current  density  is  observed  under  the  land  region  in  the 
inlet  half  where  excessive  methanol  concentration  causes  severe 
methanol  crossover  leading  to  negative  mixed  potential.  In  the 
outlet  half  region  where  methanol  concentration  is  low,  peak 
current  density  is  observed  under  the  channel  region  where 
methanol  concentration  is  high  compared  to  the  region  under 
the  land. 

3.6.  Crossover  current  density  distribution 

As  the  crossover  current  density  is  directly  determined  by 
methanol  crossover  rate  that  is  a  strong  function  of  methanol 
concentration  at  interface  between  the  anode  catalyst  layer  and 
the  proton  exchange  membrane  (see  Eq.  (15)),  calculation  result 
of  crossover  current  density  distribution  (Fig.  12)  is  very  similar 
to  methanol  concentration  distribution  (Fig.  10).  The  base  model 
has  strong  crossover  current  density  in  the  inlet  region, 
degrading  cell  performance  by  mixed  potential.  As  SGS  models 
have  reduced  current  density  in  the  inlet  region,  crossover 
current  by  EOD  is  lower  than  that  of  the  base  model  there. 
In  addition,  lowered  methanol  concentration  of  SGS  models 
also  reduces  crossover  current  by  diffusion  in  the  inlet  region. 
Flence,  SGS  models  have  reduced  overpotential  in  inlet  region. 
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Fig.  7.  Temperature  distribution  in  membrane  of  DMFCs  operating  at  i  =  150  mA  cm-2,  (a)  Base  model  (Tcen  =  313  K,  £a  =  2.4,  £c  =  2.4),  (b)  Base  model  ( Tceii  =  313  K,  £a  =  1.4, 
Cc  =  2.0),  (c)  Base  model  (Tceii  =  323  K,  £a  =  1.4,  £c  =  2.0),  (d)  SGS  model  with  variable  porosity  type  A  (Tceii  =  313  K,  £a  =  1.4,  £c  =  2.0),  (e)  SGS  model  with  tapered  channel 
(Teen  =  313  K,  Ca  =  1.4,  Cc  -  2.0). 


Average  methanol  crossover  rate  according  to  average  opera¬ 
tional  current  density  is  presented  in  Fig.  13.  SGS  models  shows 
around  3%  of  fuel  efficiency  improvement  compared  to  the  base 
model. 


3.7.  Temperature  distribution 

It  is  important  to  maintain  uniform  temperature  of  the  mem¬ 
brane  to  avoid  local  thermal  degradation,  especially  in  large  cells. 
Major  heat  sources  increasing  cell  temperature  are  current  and 
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Fig.  8.  Porosity  distribution  of  the  anode  gas  diffusion  layer  in  the  anode  channel 
direction. 


Fig.  9.  Methanol  transfer  resistance  distribution  in  the  anode  channel  direction 
(Tceii  =  313  K,  i  =  150  mA  cm"2,  ?a  =  1.4,  \c  =  2.0). 
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Fig.  10.  Methanol  concentration  distribution  in  the  anode  channel  direction 
(TCeii  =  313  K,  i  =  150  mA  cm”2,  £a  -  1.4,  £c  -  2.0). 


Fig.  12.  Crossover  current  density  distribution  in  the  anode  channel  direction 
(Tceii  -  313  K,  ~i  =  150  mA  cm'2,  £a  =  1.4,  £c  -  2.0). 


crossover  current.  Therefore,  non-uniform  distribution  of  those 
currents  leads  to  large  temperature  gradient  in  the  cell.  As  Fig.  7 
shows,  SGS  models  have  more  uniform  temperature  distribution 
than  the  base  model.  The  region  under  the  channel  at  the  inlet  is 
the  hottest  area  due  to  large  methanol  oxidation  reaction  and 
methanol  crossover. 

3.8.  Overpotential  and  cell  performance 

Overpotentials  of  each  model  are  summarized  in  Table  5.  As  all 
models  share  the  same  cathode,  cathode  overpotentials  are  almost 
the  same  for  each  case.  SGS  models  show  better  performance  than 
the  base  model  and  this  performance  improvement  came  from  the 
reduced  anode  overpotential  by  delaying  MOR  order  shift. 


Finally,  polarization  curves  are  shown  in  Fig.  14.  SGS  models 
show  improved  performance.  In  addition,  SGS  models  have  some¬ 
what  extended  limiting  current  density  because  saved  fuel  in  the 
inlet  region  due  to  reduced  MCO  was  utilized  in  the  outlet  region, 
preventing  the  cell  from  shutting  down  at  high  current  operation. 

4.  Conclusion 

It  is  found  that  DMFC  which  has  conventional  anode  structure 
suffer  from  severe  methanol  crossover  in  the  inlet  region  (low  fuel 
efficiency)  and  methanol  shortage  in  outlet  region  (low  cell 
voltage)  under  ultra-low  fuel  stoichiometry  condition.  In  addition, 
it  is  estimated  that  excessively  high  cell  temperature  degrades  cell 
performance  by  severe  methanol  crossover  in  the  inlet  region. 


Fig.  11.  Current  density  distribution  in  the  anode  channel  direction  (rcen  =  313  K,  Fig.  13.  Average  methanol  crossover  rate  according  to  average  operating  current 
i  =  150  mAcm-2,  £a  =  1.4,  £c  =  2.0).  density. 
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Table  3 

Base  cell  geometry  and  material  properties. 


Description 

Symbol 

Value 

Anode  GDL  thickness 

Adi 

270  pm 

Anode  MPL  thickness 

^mpl 

30  pm 

Anode  CL  thickness 

Acl 

62  pm 

Membrane  thickness 

^  mem 

45  pm 

Cathode  CL  thickness 

Acl 

25  pm 

Cathode  MPL  thickness 

^mpl 

30  pm 

Cathode  GDL  thickness 

Adi 

210  pm 

GDL  porosity 

£gdl 

0.6 

MPL  porosity 

£mpl 

0.4 

CL  porosity 

0.4 

GDL  permeability 

^gdl 

2.0xl0~12  m2 

MPL  permeability 

km  pi 

5.0x10  14  m2 

CL  permeability 

I<d 

1.0x10  13  m2 

Membrane  permeability 

k mem 

4.0x1 0-20  m2 

GDL  contact  angel 

Adi 

110° 

MPL  contact  angel 

$mpl 

114° 

CL  contact  angle 

Al 

96° 

Cell  length 

Aell 

75  mm 

Channel  width 

wchan 

1  mm 

Channel  height 

ffehan 

0.5  mm 

Land  width 

Wland 

1  mm 

Bipolar  plate  thickness 

Hbp 

0.5  mm 

Table  4 

Simulation  parameters. 

Description 

Symbol 

Value 

Anode  exchange  current  density 

*0,a 

45.51  A  in"2 

Cathode  exchange  current  density 

io.c 

0.12  Am"2 

Anode  catalyst  ionomer 

25  vol.  % 

Cathode  catalyst  ionomer 

25  vol.  % 

Anode  catalyst  loading 

4  mg  cm-2 

Cathode  catalyst  loading 

2  mg  cm-2 

Surface  tension 

<j 

0.0625  N  m1 

Inlet  methanol  concentration 

cMeOH 

4000  mol  itT3 

Pure  liquid  methanol  concentration 

rMeOH 

Lpure 

24719  mol  m-3 

Pure  liquid  water  concentration 

rH20 

Lpure 

55556  mol  m"3 

Total  anode  transfer  coefficients 

aa 

0.239 

Total  cathode  transfer  coefficients 

ac 

0.875 

Faraday  constant 

F 

96487  C  mol  1 

Universal  gas  constant 

R 

8.314  J  rnor’K"1 

Equivalent  weight  of  dry  membrane 

EW 

1.1  kg  mol  1 

In  order  to  overcome  these  problems  occurring  under  ultra-low 

fuel  stoichiometry  conditions, 

the  present 

study  suggested 

streamline-graded  structures  (SGS)  which 

are  designed  to  boost 

cell  performance  and  fuel  efficiency  by  mitigating  anode  non¬ 
uniformity.  Although  the  calculated  results  obtained  from  the 

present  study  showed  that  SGS  models 

achieved  about  10% 

improvement  of  voltage  and  3%  improvement  of  fuel  efficiency,  the 
better  performance  gain  can  be  obtained  through  the  optimization 

process  which  is  not  covered  in 

this  study.  In  addition,  long-term 

Table  5 

Calculation  summary. 

?a  A  Tee  11  i  (mA  cm  2)  r/a  (V) 

Vc(V) 

VCeii(V)  MCO  (%) 

Base  model  1.4  2.0  313  150 

0.444 

0.274 

0.408  12.0 

2.4  2.4  313  150 

0.396 

0.284 

0.455  19.5 

1.4  2.0  323  150 

0.477 

0.306 

0.345  16.5 

SGS  model  1.4  2.0  313  150 

0.413 

0.277 

0.443  9.7 

(Porosity  A) 

SGS  model  1.4  2.0  313  150 

0.415 

0.276 

0.441  9.9 

(Porosity  B) 

SGS  model  1.4  2.0  313  150 

0.399 

0.276 

0.458  9.0 

(Tapered) 

Average  operating  current  density  (mA  cm'2) 

Fig.  14.  Comparison  of  polarization  curves  of  DMFCs. 


stability  of  cell  operation  is  not  covered  in  the  present  study  as 
the  present  3D  model  is  a  steady-state  model.  As  this  stability  issue 
is  specially  important  for  portable  application  of  DMFCs,  it  will  be 
covered  through  the  transient  model  in  the  future. 

List  of  symbols 

Acronyms  and  abbreviations 
A  area  (m2) 

as  specific  reaction  are  (itT1) 

C  molar  concentration  (mol  m-3) 

D  diffusivity  (m2  s-1) 

Dh  Hydraulic  diameter  (m) 

E°  thermodynamic  reversible  voltage  (V) 

H  channel  height  (m) 

hm  convective  mass  transfer  coefficient 

/  current  (A) 

i  current  density  (A  m-2  or  mA  cm-2) 

i'o,a  Anode  exchange  current  density  (A  m-2) 
io,c  Cathode  exchange  current  density  (A  m-2) 

j  mass  flux  (kg  m-2  s-1),  phase  diffusion  flux 

(mol  m-2  s-1),  volumetric  transfer  current  density 
(A  m-3) 

j  phase  diffusion  flux  (mol  m-2  s-1), 

I<  permeability  (m2) 

ka  anode  kinetics  coefficient 

/<h  Henry’s  constant 

L  channel  length  (m) 

m  mass  (kg) 

M  molar  weight  (g  mol-1) 

N  molar  flux  (mol  m-2  s-1) 

p  pressure  (Pa) 

nd  electro-osmosis  drag  coefficient 

s  liquid  saturation,  entropy 

5  volumetric  source  term 

u  velocity  (m  s-1) 

U°  open  circuit  voltage  (V) 

T  temperature  (I<) 

Y  species  mass  fraction 
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Greek  letters 

a  water  transfer  coefficient,  void  fraction,  transfer  number 

e  porosity 

5  thickness 

k  ionic  conductivity 

p  density 

t  shear  stress 

y  advection  correction  factor 

r i  overpotential,  efficiency 

O  potential 

6  angle 

A  mobility,  water  content 

a  surface  tension 

Id  viscosity 

a)  catalyst  loading  amount  (g  cm-2) 


Superscript  and  subscript 

a  anode 

abp  anode  bipolar  plate 

acl  anode  catalyst  layer 

agdl  anode  gas  diffusion  layer 

ampl  anode  micro-porous  layer 

c  cathode 

capill  capillary 

cbp  cathode  bipolar  plate 

ccl  cathode  catalyst  layer 

cgdl  cathode  gas  diffusion  layer 

chan  channel 

cl  catalyst  layer 

cmpl  cathode  micro-porous  layer 

D  drag 

diff  diffusion 

e  electrolyte  phase 

eff  effective 

g  gas 


gdl  gas  diffusion  layer 
int  interfacial 

k  species 

1  liquid 

mem  membrane 

ref  reference 

s  solid  phase 

tot  total 

xover  crossover 
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